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Abstract. A resonance model for two-pion production in the pion-nucleon reaction is developed that in- 
cludes information obtained in the analysis of pion-nucleon scattering in a meson-exchange model. The 
baryonic resonances zi(1232), iV*(1440), iV*(1520), iV*(1535), and N*{1650) are included. The model re- 
produces the total cross sections up to kinetic energies of the incident pion of 350 MeV and obtains the 
shapes of the differential cross sections in reasonable agreement with the data. 

PACS. 13.75.-n Hadron-induced low- and intermediate-energy reactions - 13.75.Gx Pion-baryon interac- 
tions - 13.85.Fb Inelastic scattering:two-particle final states - 14.20.Gk Baryon resonances with S=0 



1 Introduction 

Recent experimental progress provides data for two-pion 
production in both pion-induced and electromagnetic re- 
actions up to energies of about 2.1 GeV, see e.g. Refs. PEl 
. New isobar models have been developed to deduce 
masses and widths of baryon resonances from the data [HI 
[7]. A basic problem for these analyses is given by the 
fact that the excitation of resonances is accompanied by 
other non-resonant processes, the so-called background. 
At present, the non-resonant processes are treated phe- 
nomenologically using or extending methods known 
from the analysis of pion-nucleon scattering 10, IT. 
A combined theoretical treatment of resonances and back- 
ground is a challenge for theory. 

At low energies, chiral perturbation theory provides 
a quantitative theoretical understanding of pion-induced 
two-pion production. Calculations in heavy baryon chi- 
ral perturbation theory have been extended to third or- 
der jl2lll3lll4| . An important result obtained in Ref. [T5| 
is the observation that contributions from loop diagrams 
are negligible. This is non-trivial because in pion-nucleon 
scattering, unitarity effects can be important even close 
to threshold in some partial waves. In two-pion produc- 
tion, however, the imaginary contributions interfere de- 
structively. The important contributions at third order 
are the tree level diagrams involving the finite dimen- 
sion two low-energy constants Ci and tree level correc- 
tions of order l/rriN, with mjv the nucleon mass. This 
finding of Ref. jl3) explains why a previous relativistic 
baryon chiral perturbation theory calculation at tree level 
including terms from the dimension two effective pion- 
nucleon Lagrangian works so well for pion kinetic energies 



up to 400 MeV ^Hl- Moreover, the numerical values of the 
low-energy constants (LECs) Ci of the pion-nucleon La- 
grangian can be understood in terms of resonance satura- 
tion 16j , more precisely through the s- and u-channel exci- 
tations of baryon resonances {A, A^*(1440)) and t-channel 
meson resonances {o;p). As a consequence, one may ob- 
tain a reasonable model for treating baryon resonances in 
two-pion production reactions by replacing the low-energy 
constants of chiral perturbation theory by diagrams con- 
taining explicit resonances. This has been done by Jensen 
and Miranda |^ and by Kammano and Arima ^Hl and 
in resonance models like e.g. |19II2()II21II22| . Such an ap- 
proach makes sense if one checks that at low energies the 
chiral perturbation theory results are recovered or if one 
enforces this behaviour through explicit matching of the 
pertinent amplitudes. 

In the present work, we want to include information 
about pion-nucleon scattering into the construction of a 
resonance model for two-pion production. We use results 
obtained in Ref. which treats pion-nucleon scatter- 
ing in a coupled-channel model. The model incorporates 
the effects of the ttttN states by introducing aN, ttA and 
pN channels. The interactions between these channels is 
derived from an chirally symmetric Lagrangian supple- 
mented by additional terms for the A, oj, p, and a fields. 
The model includes the Zi(1232), iV*(1520), iV*(1535), 
and iV*(1650) as explicit resonances. It is able to repro- 
duce the experimental phase shifts and inelasticities up to 
1500 MeV. For larger energies, it predicts a non- resonant 
background. The model does not need an explicit Roper 
resonance A^*(1440), but is able to describe the Pn partial 
wave and in particular the inelasticity by the dynamics of 
the aN and ttA channels. 
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Fig. 1. Classes of diagrams describing the nN — > miN reac- 
tion near threshold. The internal baryon propagators labeled 

1 and j (blue double lines) represent the nucleon, the A, and 
the N* resonances iV*(1440), iV*(1520), Af*(1535), iV*(1650), 
respectively. The p and a meson propagators are shown by 
(green) dotted lines. The (red) dashed lines stand for external 
pions, while (blue) solid lines are external nucleons. 

2 The model 

Fig- n shows the classes of diagrams which are included in 
the present resonance model. The diagrams (a), (b), and 
(c) contain two baryon propagators whereas (dl) - (d4) 
contain one. We consider the nucleon, the Z\(1232), and 
the ^*(1520), 7V*(1535), 7V*(1650), and iV*(1440) which 
are denoted by the indices i and j, respectively. We also 
include two-pion resonances in the /,r7r = and /^^r — 1 
partial waves, see diagram (d5) (and below). 

Meson-dynamical models generate unitary pion-nucleon 
scattering T-matrices by solving the Bethe-Salpeter equa- 
tion 

T ^ K + KGT, (1) 

where G denotes the two-particle propagator and K the 
scattering kernel which in principle includes all two-particle 
irreducible diagrams. The scattering kernel can be split 
into a pole term K2 which contains s-channel pole dia- 
grams and the non-pole term Ki = X — K2. As is well 
known, the T-matrix separates into a pole part T2 and a 
non-pole part Ti —T — T2 [24m25|| . The non-pole T-matrix 
is obtained by solving 

Ti = + i^iGTi . (2) 

The non-pole T-matrix generates the dressed vertex / 
which is obtained from the bare vertex /o as follows: 

/ = /o + /oGTi. (3) 

The self-energy S is given by 

S^flGf. (4) 
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Table 1. The meson-baryon interaction Lagrangian. 



vertex 




TriViV 




TvNA 


LLKAA''T^d>'nai> + h.c. 


prnv 


— gpTnTf^abcT^ad'^TTbPc 


pNN 


~9pNNi>{j^ - ^a'^''d.)Tapl^ 






aNN 




aaa 


— gacrirMcr CT^ 


TvAA 


^^A^-f5'y^Tad''7TaA'' 

-il^A^-lsYTlpl^x^ + h.c. 


pNA 


pAA 


gpNNAai-Jp O-^ d^)TaP^A 


PllTvN 


-^^^i^N-ySl^.Tad^'TTai^ + h.C. 


PllTvA 


fEll^A'^T^d^TVa^l'N' + h.C. 


PiicrA 


~gPliaNTpN'i>(^ + h.C. 




ifD^^^^^^^Yra^d,dp,'!la + h.C. 


DviTTA 


^lf^^%,-f.Tad-'TVaA^ + h.C. 


SllTvN 


^^i'M'-yf^Tad'^TVai^ + h.C. 


SiiriN 


^^^^N'-ff^d'^Vi^ + h.C. 



and the pole T-matrix reads 

T2 = Pgf, (5) 

where g^^ = G^^ — S denotes the dressed propagator. 

The diagram (dl) shown in Fig. ^ factorizes into two 
subdiagrams, the first of which is a part of the non-pole 
scattering kernel of the meson-nucleon T-matrix corre- 
sponding to a meson-exchange in the i-channel, while the 
second subdiagram is the coupling of the intermediate 
baryon propagator to a pion-nucleon final state. Likewise, 
the diagram (d2) factorizes into a pion production followed 
by a subdiagram corresponding to a non-pole meson-baryon 
scattering kernel. Diagrams (b) and (d) factorize into it- 
channel processes and pion production amplitudes. The 
M-channel and i-channel subdiagrams are the Born ap- 
proximation to the non-pole T-matrix. 

In the present work, we restrict the model to the di- 
agrams shown in Fig. ^ This allows a relatively simple 
treatment of two-pion production because only tree-level 
diagrams with dressed vertices and dressed propagators 
taken from a pion-nucleon scattering model have to be 
evaluated. An iteration of the non-pole subdiagrams would 
require to solve Eq. Q which commonly is done numeri- 
cally relying on a partial wave representation in the center- 
of-mass frame, but implies technical complications when 
boosts to other frames are required, e.g. when evaluating 
the two-pion production amplitude. 

In Ref. I^ni, the inelasticities of the Pn partial wave 
were explained by the final-state interactions, i.e. the non- 
pole contribution of the T-matrix. In the present model, 
such a structure cannot be generated. We therefore have 
to treat the Roper in a simplified way and represent it as 
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Table 2. Coupling constants and masses in GeV. The con- 
stants gi and g2 are in GeV~^. 



/^iviv/(47r) 


0.0778 


./pil{1440)7r]v/(4'?'") 


0.011 


/^]vzi/(47r) 


0.36 


./pil{1440)7rZi/(4'''") 


0.04 


Ppiv jv/(47r) 


0.80 


/pil{1440)<Tjv/(47i") 


13.0 




1.94 


./D13(1520)7riv/(4'r) 


0.0009 




1.03 


/oi3(1520)7rZi/(4'I") 


0.03 




0.04 


/sil{1535)7r]v/(4''I") 


0.003 


fpNA/i^^) 


4.5 


/si1{1535)t,]v/(4'''") 


0.47 


Ppzlzi/(47r) 


16.0 


/sil(1650)r,]v/(47I") 


0.009 


Hp A A 


15.0 






gpTT^/ (47r) 


2.905 


ff?/(47r) 


98.94 


5L<./(47r) 


0.625 


52V(47r) 


7.32 


mjv 


0.93893 


"T'I313{1520) 


1.515 




1.232 


"1-511(1535) 


1.535 


'^Pll{1440) 


1.491 


"1-311(1650) 


1.701 




0.13803 


Mp 


0.772 




0.5473 




0.8346 



an s-channcl resonance. A Roper propagator is introduced 
which has a structure analogous to the nucleon propaga- 
tor. 

We simplify the formalism employed in Ref. by 
solving the coupled channel problem for meson-nucleon 
scattering in the X-matrix approximation. Moreover, we 
employ a derivative coupling for the coupling of the nu- 
cleon resonances to the pseudo-scalar mesons which was 
used in Ref. [211 to improve the description of the SuttN 
partial wave. 

For brevity, we only display the various propagators 
that enter the calculation. All the necessary formalism to 
calculate the relevant amplitudes, cross sections etc. can 
e.g. taken from Ref. 

The nucleon propagator is given by: 



Sn{p) 



TUN 



(6) 



with niN the nucleon mass. The self energy Sn is obtained 
by including the ttN, ttA, and the aN reaction channels 
in Eq. ||2J). The propagator of the zi(1232) is given by 



{P 



(7) 



with the mass of the Z\(1232). The spin-1/2 contribu- 
tion is identical with the spin-1/2 part of the Rarita- 

Schwinger propagator. These contributions are, however, 
non-propagating and can be represented in an effective 
field theory approach by contact operators, see e.g. |28j . 
The delta self-energy is evaluated taking into account 
only the pion-nucleon intermediate state and neglecting 
the non-pole contribution of the pion-nucleon interaction, 
i.e.: 

Sa - fl G^N h ■ (8) 



This is a good approximation because the P33 partial wave 
in the pion-nucleon scattering model of Ref. 23 is given 
mainly by the Z\-pole diagram for partial waves up to 
approximately 1.3 GeV. 

The interaction Lagrangians employed are shown in 
Tabled A partial refit of the parameters of Ref. [221 is re- 
quired. In a first step, the parameters of the TriV potential 
are readjusted without the N* pole diagrams to the scat- 
tering lengths and phase shifts below 1.2 GeV. Then the 
ttAA, pNA, and the pAA couplings are fitted. Finally the 
parameters of the N* resonances are fixed. Table |21 sum- 
marizes the coupling constants and masses. The value of 
the tensorial pNN vertex has been reduced in comparison 
to the value employed in Ref. 

The diagrams (dl), (d2), (d3), (d4), and (d5) con- 
tain a cr or p propagator. Since here we need a reason- 
able model for the pion-pion interaction below two-pion 
invariant masses of 1 GeV only, a simplified version of 
the meson-exchange model to pion-pion scattering suf- 
fices [27|- We include a dressing of the a and p propaga- 
tors by two-pion intermediate states only, neglecting KK 
intermediate states and non-pole scattering kernels. Ta- 
bles ^ and |21 show the relevant interaction Lagrangian and 
coupling constants. Unitarity is taken care of by solving 
Eq. (QJ for the two-pion system. 



3 Results 

In Fig. 13 the threshold behaviour of the total cross sec- 
tions of the five different experimentally accessible chan- 
nels 



TT^ p 
TT^ p 
TT^ p 



TT^ Tt" P , 



(9) 



are shown as a function of the kinetic energy of the initial 
state pion in the laboratory frame, T^-. The corresponding 
invariant mass ^/s of the initial pion-nucleon system is 
given by 



{rriN + Mjr) + 2mArr^ 



(10) 



with the relevant neutral or charged pion mass de- 
pending on the channel under consideration. The Z\-isobar 
corresponds to T^r = 0.19 GeV, while the maximal pion ki- 
netic energy = 0.4 GeV translates into y/s = 1.38 GeV 
which is below the nominal mass of the Roper resonance. 
The observables are calculated assuming isospin symme- 
try. Following Refs. [OIIT^ , we included the effect of isospin 
breaking by the masses of the final states by shifting the 
isospin symmetric threshold to the correct threshold en- 
ergy for each reaction channel. 

Close to threshold, the present model obtains cross sec- 
tions which agree with the results of chiral perturbation 
theory. This could be expected since the model incorpo- 
rates the experimental information on low energy pion- 
nucleon scattering. It should be stressed, however, that 
the model is not as precise as chiral perturbation theory 
because it does not offer a counting scheme which would 
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0.2 0.3 0.4 



TJGeV] 

Fig. 2. Total cross sections for the reaction nN nnN. The 
(red) solid lines show the results of the present model. For 
comparison, the results obtained in Heavy Baryon chiral per- 
turbation theory |13| are displayed by the (green) dot-dashed 
line. The data are taken from Ref. and the compilation in 

Ref. unj. 



allow systematic improvements. Chiral perturbation the- 
ory starts to overshoot the experimental cross sections for 
the TT+p TT+vr+n reaction at about ~ 0.3 GeV. In- 
creasing the order of the expansion would push the range 
of validity of the chiral calculation to higher energies. In 
the model, some higher order terms are included by uni- 
tarization effects. It is important to note that these higher 
order effects are of no relevance in the threshold region. 
The model continues to agree with the experimental data 
up to about T„ ~ 0.4 GeV. 

In Fig. O we investigate the effect of different reso- 
nances on the total cross sections by omitting the con- 
tribution of various baryon resonances. The Z\- resonance 
may be expected to play an important role for the Tr~^p 
reaction for pion kinetic energies of about = 0.2 GeV. 
This is seen in the reaction Tr~^p —> tt'^tt'^p, but not in 
the reaction ir^p — )■ 7r+7r+?i. A closer inspection showed 
that the contribution of a diagram with an intermediate 




0.2 0.3 0.4 



TJGeV] 

Fig. 3. Influence of resonances on the total cross sections for 
the reaction nN mrN. Full model: (red) solid lines; omis- 
sion of the A resonance: (green) dot-dashed lines; omission of 
the A'^* resonances: (blue) dashed lines. The (violet) dotted 
lines show the contribution of diagrams with an intermediate 
p-meson decaying to two pious. 

A resonance by itself is large, but that there is destructive 
interference between the A contributions. Such cancella- 
tions are also observed in the description of electromag- 
netic two-pion production off nucleons in the threshold 
region, see e.g. [21]. The effect of the A can be seen in the 
TT~p reaction channels. The impact of the N* resonances 
5ii(1535) and _Di3(1520) on the total cross sections is neg- 
ligible in each of the five reaction channels in the energy 
region investigated. The largest effect of N* resonances 
is due to the Roper resonance and can be seen in the 
Tr~p — > TT+TT^n and Tr~p — > Tr'^Tr^n reactions, where the 
two final pions can be produced in a relative s-wave. The 
role of the Roper was recently also investigated in Ref. |31j . 
The aN channel is not allowed in the n^p — > t:^t:~p re- 
action and therefore the TV* resonances contribute very 
little. 

The model includes two-pion production via an in- 
termediate rho-meson. For the 7r+7r+n and Tr^Tr^n final 
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Fig. 4. Differential cross sections dcr/dM^^ for TT+p —> ■K'^i^^n. Fig. 5. Differential cross sections da/dt for -K^p — > -K^Ti^n. 

Full model: (red) solid line; phase space: (blue) dotted line; Full model: (red) solid line; phase space: (blue) dotted line; 

chiral perturbation theory: (green) dot-dashed line 1131 . The chiral perturbation theory: (green) dot-dashed line |13| . The 

data are from Ref. 0- data are from Ref. In- 



states, trivial isospin selection rules do not allow contri- 
butions from an intermediate rho-meson, and in the other 
final states, the effect of those processes is marginal, as 
expected, see Fig. |31 

The differential cross sections da/dM^^ for the reac- 
tions TT^p —^ TT^TT^n are shown for five different values of 
the kinetic energy T.^ of the incident pion in Fig. ^ (here 
Mtttt is the invariant mass of the final-state two-pion sys- 
tem). The data measured at TRIUMF can be fitted rea- 
sonably well by three-body phase space [5]. Chiral per- 
turbation theory (CHPT) predicts the average magnitude 
of the cross sections correctly, but produces a large devi- 
ation from phase space emphasizing low values of M^^. 
The cross sections calculated within the resonance model 
agree with chiral perturbation theory for low values of 
for both T^^ = 223 MeV and T^^ = 243 MeV, but 
drop with increasing M^^ more slowly than the CHPT 
predictions do. For — 305 MeV, the model reproduces 
the experimental cross sections above M^^ = 6M^, but 



overestimates the data below M^^ = 6M^. We note again 
that the higher order terms generated by the unitarization 
procedure play only a small role at low energies. 

In contrast to the M-^tt distributions shown in Fig. 0] 
the differential cross sections da/dt strongly deviate from 
phase space, see Fig. |31 Chiral perturbation theory repro- 
duces the shapes of the ^-distributions rather well, but 
misses the strong rise of the data at small momentum 
transfers. The model can reproduce the t-distributions 
nearly quantitatively for small values of the momentum 
transfer t, but overpredicts the data at large momentum 
transfers. 

The differential cross section da /d cos 6 for the 7r+p — > 
TT+TT+n reaction has to be symmetric for cos 9 = because 
of the symmetry of the reaction under exchange of the two 
produced pions. The experimental cross sections show a 
maximum at cos 9 = 0, see Fig. The chiral perturba- 
tion theory at order three predicts a minimum, however. 
The model obtains a maximum near 9 = 90° in qualita- 
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Fig. 6. Differential cross sections da/dcos9 for Tr^p 
TT+TT^n. Full model: (red) solid line. A calculation employing 
bare vertices /o is displayed by the (blue) dashed line. The pre- 
dictions of chiral perturbation theory are shown by the (green) 
dot-dashed line |13|. The data are from 0. 
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Fig. 7. Differential cross sections da/dM^^ for Tv~p Tr'^Tv~n. 
Model without Roper resonance: (red) solid line; Model includ- 
ing Breit-Wigner Roper resonance: (red) double-dotted dashed 
line; phase space: (blue) dotted line; chiral perturbation the- 
ory: (green) dot-dashed line JHl- The data are from Ref. 0- 




tive agreement with the data. When replacing the dressed 
vertex functions / by the bare ones, however, the shape of 
the angular distribution changes and one obtains a mini- 
mum at 9 = 90°. This finding suggests that a chiral per- 
turbation theory carried to fourth order would solve the 
problem of the angular distributions in the n^p 7r+7r+n 
reaction. 

The differential cross sections da/dM^^ for the reac- 
tion for TT^p TT+TT^n obtained at TRIUMF are shown in 
Fig. 13 The experimental data deviate significantly from 
phase space and show a pronounced maximum close to 
the largest values of M^^ . Chiral perturbation theory re- 
produces the asymmetric distributions very well for the 
pion kinetic energies = 223 MeV and = 243 MeV. 
The model reproduces the shapes of the distributions, but 
overestimates the absolute magnitudes. When switching 
off the contribution from the Roper resonance, however, 
a reasonable reproduction of the data is achieved. In the 



present model, we have treated the Roper in a simpli- 
fied way, which apparently produces an overestimation of 
the effect of that resonance. If one wants to incorporate 
the possibility to generate resonances dynamically, a ma- 
jor revision of the present approach is required. In order 
to include the full non-pole T-matrix, an iteration of the 
non-pole scattering kernel in the two-pion production di- 
agrams has to be performed for different Lorentz frames. 
This goes beyond the scope of the tree-level like model 
discussed here. 



4 Summary 

We have presented a resonance model for the ttN — + ttttA^ 
reaction which incorporates information from pion-nucleon 
scattering - these two processes are intimately connected 
and should not be treated independently from each other. 
The model is able to reproduce the total cross sections 
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for kinetic energies of the incident pions up to about 350 
MeV for all five reaction channels. This is a success of the 
present approach which is only partially shared in spe- 
cific reaction channels by other resonance models which 
do not include consistently pion-nucleon scattering and 
two-pion production. The agreement with the differential 
experimental cross sections is good, though not perfect. 
The total cross sections have a smooth dependence on the 
pion kinetic energy which does not allow to deduce ob- 
vious information about resonances. The differential cross 
sections, on the other hand, show characteristic deviations 
from three-body phase space. The inclusion of the A res- 
onance suffices to explain the overall shapes. 

In the second resonance region, the model starts to 
break down. Here, we could trace the difficulties to our 
approximate treatment of the Roper resonance. The re- 
sults obtained in the present study suggest to go beyond 
the Born approximation of the non-pole pion-nucleon T- 
matrices, when constructing the two-pion production am- 
plitudes for analyses in the second resonance region. 



Acknowledgements 

We thank Christoph Hanhart for discussions. Partial sup- 
port from the EU Integrated Infrastructure Initiative Had- 
ron Physics Project (contract no. RII3-CT-2004-506078) 
and DFG (SFB/TR 16, "Subnuclear Structure of Mat- 
ter" ) is gratefully acknowledged. 



References 

1. R. Fatemi et al. [CLAS Collaboration], Phys. Rev. Lett. 91 
(2003) 222002 arXiv:nucl-ex/0306019 . 

2. W. Langgartner et al, Phys. Rev. Lett. 87 (2001) 052001. 

3. M. Wolf et al, Eur. Phys. J. A 9 (2000) 5. 

4. S. Prakhov et al. [Crystal Ball Collaboration], Phys. Rev. 
C 69 (2004) 045202. 

5. M. Kermani et al. [CHAOS Collaboration], Phys. Rev. C 
58 (1998) 3419. 

6. I. G. Aznauryan, V. D. Burkert, G. V. Fedotov, 
B. S. Ishkhanov and V. I. Mokeev, Phys. Rev. C 72 (2005) 
045201 arXiv:hep-ph/0508057 . 

7. A. Anisovich, E. Klempt, A. Sarantsev and U. Thoma, Eur. 
Phys. J. A 24 (2005) 111 arXiv:hep-ph/0407211 . 

8. T. P. Vrana, S. A. Dytman and T. S. H. Lee, Phys. Rept. 
328 (2000) 181 arXiv:nucl-th/9910012 . 

9. R. A. Arndt, W. J. Briscoe, I. I. Strakovsky and 
R. L. Workman, Phys. Rev. C 66 (2002) 055213 
farXiv:nuc l-th/0205067 . 

10. L. Tiator, D. Drechsel, S. Kamalov, M. M. Giannini, 
E. Santopinto and A. Vassallo, Eur. Phys. J. A 19 (2004) 
55 arXivmucl-th/0310041 . 

11. T. Feuster and U. Mosel, Phys. Rev. C 58 (1998) 457 
farXiv:nucl-th/9708051 . 

12. V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. B 
457 (1995) 147 arXiv:hep-ph/9507418 . 

13. N. Fettes, V. Bernard and U.-G. Meifiner, Nucl. Phys. A 
669 (2000) 269 arXiv:hep-ph/9907276| . 



14. N. Mobed, J. Zhang and D. Singh, Phys. Rev. C 72 (2005) 
045204. 

15. V. Bernard, N. Kaiser and U.-G. Meil3ner, Nucl. Phys. A 
619 (1997) 261 arXiv:hep-ph/9703218 . 

16. V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. A 
615 (1997) 483 arXiv:hep-ph/9611253 . 

17. T. S. Jensen and A. F. Miranda, Phys. Rev. C 55 (1997) 
1039. 

18. H. Kammano and M. Arima, Phys. Rev. C 69 (2004) 
025206. 

19. E. Oset and M. J. Vicente- Vacas, Nucl. Phys. A 446 (1985) 
584. 

20. O. Jaekel, M. Dillig and C. A. Z. Vasconcellos, Nucl. Phys. 
A 541 (1992) 675. 

21. J. A. Gomez Tejedor and E. Oset, Nucl. Phys. A 600 
(1996) 413 arXiv:hep-ph/9506209 . 

22. J. C. Nacher, E. Oset, M. J. Vicente and L. Roca, Nucl. 
Phys. A 695 (2001) 295 arXiv:nucl-th/0012065 . 

23. O. Krehl, C. Hanhart, S. Krewald and J. Speth, Phys. Rev. 
C 62 (2000) 025207 arXiv:nucl-th/9911080 . 

24. R. W. Haymaker, Phys. Rev. 181 (1969) 2040. 

25. I. R. Afnan and A. T. Stelbovics, Phys. Rev. C 23 (1981) 
1384. 

26. A. M. Gasparyan, J. Haidenbauer, C. Hanhart 
and J. Speth, Phys. Rev. C 68 (2003) 045207 
arXivmucl-th/0307072 . 

27. F. P. Sassen, S. Krewald and J. Speth, Phys. Rev. D 68 
(2003) 036003 arXiv:hep-ph/0212056 . 

28. V. Bernard, T. R. Hemmert and U.-G. Meifiner, Phys. 
Lett. B 565 (2003) 137 arXiv:hep-ph/0303198 . 

29. V. V. Vereshagin et al., Nucl. Phys. A 592 (1995) 413 
arXiv:hep-ph/9504361 . 

30. V. Bernard, N. Kaiser, U.-G . Mei finer and A. Schmidt, 
Nucl. Phys. A 580 (1994) 475 "arXrv:nucl-th/9403013j. 

31. H. Kamano and M. Arima, arXivmucl-th/0601057r "^ 



